METHOD AND DEVICE FOR RECEIVING A RADIO SIGNAL 



BACKGROUND OF THE INVENTION 

The present Invention relates to digital radiocommunication techniques 
using Code Division Multiple Access (CDMA). 

A CDMA signal processed by a receiver can be expressed, after 
filtering and transposition to baseband, as follows: 

y(t)=i;y (1) 

u=1 

where n(t) is additive noise and U is the number of channels multiplexed on the 
CDMA carrier, whose contributions y"(t) have as a general expression: 

* y^t) = Xbrs[^(tHT) (2) 

where: 

• bH is the information symbol of rank i transmitted on the u-th channel; 

• sH(t) is a generalized code given by the convolution of the impulse 
response of the u-th channel with the portion corresponding to the symbol 
of the spreading code c" assigned to the channel. 

The number U corresponds to the number of user if each user involved 
utilizes a single channel. There may however be several channels per user. 

The spreading codes c" are sequences of discrete samples called 
"chips", with real values (±1) or complex values (±1±j), having a given chip rate. 

The symbols b" also have real values (±1) or complex values (±1±j). The 

duration of a symbol on a channel is a multiple of the chip duration, the ratio of 
the two being the channel spreading factor Q. 

In certain systems, the spreading factor may vary from one chainnel to 



another. In such a case, a common spreading factor Q is considered, equal to 
the greatest common divisor (GCD) of the U spreading factors Q". A symbol on 
the channel u is then regarded as the concatenation of Q"/Q consecutive 
symbols bH whose values are identical. 

The duration of the generalized response sH(t) corresponds to Q+W-1 

chips if W denotes the length of the impulse response expressed as a number 
of chips. 

By sampling at the chip rate the CDMA signal y(t) received for a block 
of n symbols on each of the channels, the receiver obtains complex samples 
that can be modeled by a vector Y of nxQ+W-1 components: 

Y = A.b + N (3) 

where: 

• b denotes a column vector of size nxU, which can be decomposed into 

=(b![,b2, ",bj), where represents the transposition operation, the 
vectors bj being of size U for 1 ^i ^n, with bT = (b?,bp,- ^bf^); 

• N is a random noise vector of size nxQ+W-1; 

• A = (A.|, Ag, .... Ap) is a matrix of generalized codes of size (nxQ+W-1)x 
(nxU) which can be subdivided into n sub-matrices Aj of size (nxQ+W-1) 
xU. In the matrix Aj (1 ^i ^n), the u-th column (1 is a convolution of 
the impulse response of the u-th channel and of the Q samples of the 
spreading code of the u-th channel corresponding to the i^th symbol of 
the block. 

In other words, the matrices Aj may be written: 

Ai=(Qj'.n2,....Qu) (4) 

with: QH=mH.HH (5) 
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where is a Toeplitz matrix of size (nxQ+W-1) x (n>«Q+W-Q) obtained from 

the values cH(q) of the chips of the spreading code c" of the u-th channel over 
the duration of the i-th bit of the block: 



M,H = 



cH(1) 0 
cH(2) Ci"(1) 
Ci"(2) 

cH(Q) 
0 cf(Q) 



0 

c[^(2) 
0 c"(Q) 



(6) 



and hH is a column vector of size (n-1)xQ+w which, when the U channels are 
received synchronously, contains (i-1)xQ zero components, followed by the W 
samples of the u-th channel's Impulse response relating to the i-th symbol bH , 
and followed by (n-i)xQ other zero components. The time offsets in reception 
along the various channels, In numbers of chips, are manifested as 
corresponding offsets of the W samples of the impulse response of the 
channels along the vector HH . 

The receiver most commonly used, the so-called rake receiver, uses 
one or more matched filters to estimate the value of the symbols transmitted on 
each channel from an estimate of the impulse response of the channel along 
one or more propagation paths. 

The operation performed by such receivers amounts to performing the 
matrix product: 



Z = A*.Y 



(7) 



where A* is the conjugate transpose of an estimate A = (A^, A2 A^) of the 

matrix A, the matrices Aj stemming from the impulse responses estimated by 
applying relations (4) and (5). 
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The nxu components of the vector Z are respective soft estimates 
of the nxu symbols bj" of the vector b. If the decodings performed downstream 
admit soft estimates as input, the components of the vector Z can be used 
directly. Otherwise, the sign of these components is taken to form the hard 
5 estimates of the symbols. 

The matched filter receiver is optimal when the generalized codes 
(vectors Q.") are pairvyise orthogonal, i.e. when the correlation matrix R = A*.A 
is diagonal. In general, the systems adopt pairwise orthogonal spreading codes 
having good autocorrelation properties, whereby this condition is fulfilled to a 
10 first approximation. 

However, when the impulse response of the channel is taken Into 
account, the orthogonality condition is no longer fulfilled. The above 
approximation becomes poor especially in the presence of multiple propagation 
paths. 

15 Certain receivers cany out a posteriori correction of the soft estimates 

of the symbols emanating from the matched filter receiver by taking account of 
the inter-users interference and/or of the inter- symbol interference on one and 
the same channel, thereby substantially improving the performance. Such a 
procedure, based on a so-called MFPIC ("Matched Filter Parallel Interference 

20 Cancellation") algorithm, is disclosed in WO 01/99301. Its advantage is that it 
does not overly increase the global complexity of the calculations with respect 
to the traditional "rake" receiver. However, these a posteriori corrections follow 
upon an optimization of the system (3) relying on the above approximation. 

This MFPIC algorithm belongs to the class of multiuser detection 
!5 algorithms (MUD) that offer better perfonnance than the regular matched filter 
receiver. This class includes more accurate algorithms than MFPIC, in that they 
can take into account the terms of the matrix R that are farther from its 
diagonal. 

A good example of a MUD algorithm useable in such a context is the 
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SDP algorithm described in the article by M. Abdi, et al., "Semidefinite Positive 
Relaxation of the Maximum-Likelihood Criterion Applied to Multiuser Detection 
in a CDMA Context". IEEE Signal Processing Letters, Vol. 9, No. 6, June 2002, 
pp. 165 167. 

5 The better performance of these MUD algorithms is obtained at the 

price of greatly increased complexity. In general, their complexity is more than 
linear in the size nxU of the problem to be solved, this rendering them very 
expensive by comparison with more conventional receivers such as the "rake" 
or the MFPIC. 

10 An object of the present invention is to find a good compromise 

between the performance and the complexity of a CDMA receiver. 

SUMMARY OF THE INVENTION 

The invention thus proposes a method of processing a signal received 
via a radio interface and including contributions from a plurality of channels 
15 multiplexed by respective spreading codes. The method comprises the steps 
of: 

/a/ estimating response parameters of the multiplexed channels; 
Ibl calculating soft estimates of symbols transmitted over the multiplexed 
channels, as a function of the received signal and of the estimated 
20 response parameters; 

Id dividing the symbols whose soft estimates have just been calculated 
between a first set of symbols satisfying a confidence criterion applied to 
said soft estimates and a second set of symbols not satisfying the 
confidence criterion; 

25 /d/ determining a modified signal by subtracting estimated contributions 
corresponding to the symbols of the first set, respectively, from the signal 
subjected to the previous calculation of soft estimates; and 
lei calculating new soft estimates of the symbols of the second set only, as a 
function of the modified signal and of the estimated response parameters. 
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The symbols for which the first calculation affords sufficiently reliable 
estimates are no longer to be estimated in the second calculation. Their 
contribution to the signal is deducted so as to take into account the interference 
that they generate for the other symbols. 

5 Various detection algorithms can be employed for the successive 

calculations of soft estimates. 

In a preferred embodiment of the method, step lei is executed 
according to a detection algorithm of more complex nature than step /b/, in 
particular according to a multi-user detection algorithm. The method then profits 
10 from the performance of a complex algorithm of MUD type, but by applying the 
latter to a problem of smaller size than the problem forming the subject of the 
first calculation. 

The method thus makes it possible to achieve a compromise between 
the symbols detection performance and the complexity of the required 
15 calculations. In the design or the configuration of the receiver, it is possible to 
favor the performance or the complexity in this compromise, by appropriate 
adjustment of the confidence criterion employed. 

It is moreover possible to repeat the sequence of steps lol to lei one or 
more times. The larger or smaller number of iterations and/or the severity of the 
20 confidence criterion applied at each iteration (this confidence criterion may vary 
from one iteration to the next) are parameters that may also be adjusted as a 
function of the performance/complexity compromise sought. 

Another aspect of the present invention pertains to a device for 
processing a signal received via a radio interface, including contributions from a 

25 plurality of channels multiplexed by respective spreading codes. This device 
comprises means for estimating response parameters of the multiplexed 
channels, first means for calculating soft estimates of symbols transmitted over 
the multiplexed channels, as a function of the received signal and of the 
estimated response parameters, means for dividing the symbols between a first 

30 set of symbols satisfying a confidence criterion applied to the soft estimates 
produced by the first means for calculating and a second set of symbols not 
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satisfying the confidence criterion, means for determining a modified signal by 
subtracting estimated contributions corresponding to tlie symbols of the first 
set, respectively, from the received signal, and second means for calculating 
new soft estimates of the symbols of the second set only, as a function of the 
5 modified signal and of the estimated response parameters. 

Such a device can in particular be incorporated Into a base station of a 
CDMA radiocommunication system. 

When several iterations are envisaged, the device further comprises 
second means for dividing the symbols of the second set between a first subset 

10 of symbols satisfying a second confidence criterion applied to the soft 
estimates produced by the second means for calculating and a second subset 
of symbols not satisfying the second confidence criterion, means for 
determining a second modified signal by subtracting estimated contributions 
corresponding to the symbols of the first subset, respectively, from the modified 

15 signal which was subjected to the second means for calculating, and third 
means for calculating new soft estimates of the symbols of the second subset 
only, as a function of the second modified signal and of the estimated response 
parameters. 

Another aspect of the present invention pertains to a computer program 
20 to be installed in a radiocommunication receiver, the program comprising 
instructions for implementing a method as defined hereinabove during 
execution of the program by a signal processing unit of the receiver. 

BRIEF DESCRIPTION THE DRAWINGS 

Figure 1 Is a schematic diagram of an exemplary reception device 
25 according to the invention. 

Figure 2 is a flow chart of an exemplary method according to the 
invention. 

DESCRIPTION OF PREFERRED EMBODIMENTS 



The device represented in Figure 1 is part of tlie receiving stage of a 
radiocommunication station able to communicate witli several remote stations 
1. 

The uplink channels used by these distant stations 1 are multiplexed by 
the CDMA technique, so that the radio signal picked up by the antenna 2, and 
converted to baseband, can be represented in the form (1)-(2) for U multiplexed 
channels originating from V stations (1 <V ^U). 

The station incorporating the device is for example a base station of a 
third-generation cellular radiocommunication system of the UMTS ("Universal 
Mobile Telecommunication System") type. 

In Figure 1. the unit 3 diagrammatically represents the modules 
performing in a conventional manner the signal reception preprocessing 
(amplification; filtering, conversion to baseband, sampling at the chip 
frequency). This unit 3 delivers blocks Y of nxQ+W-1 samples, corresponding 
to blocks of n symbols transmitted simultaneously on the U channels. If the 
blocks of n symbols follow one another without interruption on the channels, 
there is an overlap of W samples (chips) between the successive blocks Y. 
corresponding to the duration of the impulse response. 

The received-signal blocks Y are provided to a module 4 which 
estimates the impulse responses IHH of the U multiplexed CDMA channels, 
with the aid of correlations with the spreading codes produced by 
pseudorandom code generators 5 (1 ^u^U. 1 <i<n). The module 6 then 
estimates the matrix A of the generalized codes according to (4) and (5), i.e. its 
[(i-1)xU + u]-th column is given by qH =MH.HH. 

With the aid of the parameters of the matrix A, a first calculation of soft 
estimates is applied to the block Y. In the example represented, this first 
calculation is effected by the modules 7 to 9 according to the MFPIC algorithm 
described in WO 01/99301. The module 7 performs a conventional detection of 
"rake" type on each channel, according to relation (7) above. It produces first 
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soft estimates ZH of the transmitted symbols b" . These first estimates zH 
could be used directly, but their representativity is improved by refining them in 
the module 8, which applies the correction envisaged In the MFPIC algorithm to 
take account of inter-symbol interference and/or inter-user interference. In the 
5 case where the symbols transmitted are bits, this correction is expressed by: 

sf,"=2f-|:Ri;j".sgn(z;')-i:tRjj^sgn(z)') (8) 

j=1 V=1j=1 

Where sgn(.) designates the sign function with values in {-1,+1} and RHj^ 
designates the term situated in the [(1-1 )xu + u]-th row and the [0-1 )xU + v]-th 
column of the correlation matrix R = A*.A calculated by the module 9. 

"•o The first term subtracted in (8) corresponds to inter-symbol interference 

on the u-th channel while the second term subtracted corresponds to the inter- 
channel Interference. The estimates sf," produced by the module 8 are 

approximations of the "softblts", to within a multiplicative coefficient 4/o^, where 
a designates the power of the additive noise picked up, an estimate of which is 
15 conventionally made available by the probing module 4. The sign of sfj" forms 

a hard estimate of the bit bH , while its absolute value measures the likelihood 
of this estimate. 

According to the invention, these soft estimates sf " are examined by a 

module 10 In such a way as to identify a set F of symbols whose estimates are 
20 regarded as the most reliable. 

By way of example, the module 10 sorts the estimates sfj" in the order 
of descending absolute values and places in the set F the K% of the symbols 
whose estimates are ranked first. Another possibility Is to place In the set F the 
symbols whose estimates sf" have an absolute value greater than a 
25 confidence threshold p that can be taken proportional to 4/<y^. The threshold p 
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or the percentage K can be made adaptive, for example as a function of the 
signal-to-noise ratio observed by the receiver. 

For the symbols of the set F which is thus determined, the device will 

deliver the soft estimates sf " produced by the module 8 or hard estimates b" 

consisting of their signs, according to the requirements of the processing 
circuits situated downstream. 

The symbols that do not form part of this set will form the subject of a 
second calculation of soft estimates which preferably uses an MUD algorithm. 

The columns of the matrix of the generalized codes A are sorted by a 
module 1 1 . Each column associated with a symbol of the set F is extracted 
from A so as to be multiplied by the hard estimate bH of this symbol (multiplier 
12), thereby giving an estimate of the contribution of the interference caused by 
this symbol on the others. This contribution bH.QH is subtracted from the vector 

Y representing the input signal by the subtracter 13. Such a subtraction is 
performed for each symbol of the set F: 

Y' = Y- Z^i -^f (9) 
(i.u)eF 

= Zi;b[^^"+N- ^S^nf (10) 
i=1u=1 (i.u)eF 

« Zb-'.^i'+N (11) 
(i.u)eF 

The relation (11) shows that by taking into account the contributions of 
the properly estimated symbols, the size of the system to be processed is 
reduced from nxU to nxU - card(F). The size reduction Is very appreciable if it 
has been possible to properly estimate sufficient symbols during the first pass. 
This facilitates recourse to an algorithm of more complex nature for the second 
estimate of the remaining symbols, in particular to an MUD algorithm. 
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This second estimate is performed by the detection module 14 of 

Figure 1, to which the modified signal V is subjected. Each column vector qH 

with (i,u) G F is removed from the matrix A by the module 11 so as to form a 
reduced matrix of generalized codes according to which the module 15 obtains 

5 the reduced correlation matrix R = A*.A. The vector Y' and the reduced 
matrices A and R constitute the input data for the MUD algorithm implemented 
by the module 14, which is for example the aforesaid SDP algorithm. 

The soft estimates sf." (or hard estimates) produced by the MUD 

module 14 are ultimately combined with the estimates determined in a 
10 sufficiently reliable manner by the MFPIC algorithm in the course of the first 
pass, so as to reconstruct the set of demodulated data. 

The soft estimates sf|" produced by this MUD module 14 may also 

form the subject of a sorting between reliable estimates and unreliable 
estimates, the latter then forming the subject of a new calculation of estimates. 
15 This sorting and re-estlmation process can be repeated a certain number of 
times. Such iterative procedure is illustrated by Figure 2. 

The first two steps 20 and 21, prior to the iterations, consist of the 
conventional estimation of the responses of the channels and of the 

generalized codes qH (operations of the modules 4 and 6 of Figure 1). In step 

20 22, the set E of the symbols to be estimated in the course of the next iteration 

is initialized to the Cartesian product {1,2 n} x {1,2,...,U} corresponding to 

the entirety of the symbols transmitted within the current block over the 
channels to be processed. 

In each iteration, the first step 23 consists In assembling the column 
25 vectors qH associated with the symbols of the set E to form the matrix of 
codes A which will be used for the next detection, and in calculating the 
correlation matrix R = A*.A . The detection is effected In the next step 24 on 
the basis of the block Y with the aid of the matrices A and R . By way of 
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example, the algorithm employed In step 24 is the MFPIC in the first Iteration 
and the SDP in each subsequent iteration. 

The soft estimates obtained In step 24 are subjected to the confidence 
criterion so as to construct, in step 25, the set F of indices (i,u) of the most 
. reliable symbols, which are for example the Indices (l,u) of the set E such that 

|sfj"| ^p. 

If all the estimates are deemed sufficiently reliable (F = E In test 26). 

the detection procedure Is terminated and the soft estimates sfj" which were 

calculated may be delivered In step 27. The threshold p can vary In the course 
of the iterations. It may in particular decrease so as to make the confidence 
criterion less and less severe. To limit the procedure to a maximum number X 
of iterations, it Is possible to set p = 0 for the X-th iteration. 

If certain estimates are not sufficiently reliable {F E \r\ test 26), the 
hard estimates b" of the symbols of F are determined In step 28. The 
estimated contribufions bH.QH of these symbols are subtracted from block Y in 
step 29 (operafion of the modules 12 and 13 of Figure 1). To Initialize the next 
iteration, the set E of symbols to be estimated Is updated In step 30 by deleting 
therefrom the correctly estimated symbols of the set F. 

A procedure such as that Illustrated by Figure 2 can be implemented by 
programming a digital signal processor provided In the radio receiver. 



